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Abstract 

We investigate direct detection rates for Dark Matter candidates arise in a SU(2)l x U(1)y 
with an additional doublet Higgs proposed by Barbieri, Hall and Rychkov. We refer this 
model as 'Heavy Higgs Model'. The Standard Model Higgs mass comes out in this model very 
heavy adopting the few per cent chance that there is no Higgs boson mass below 200 GeV. 
The additional Higgs boson develops neither any VEV due to the choice of coefficient of the 
scalar potential of the model nor it has any coupling with fermions due to the incorporation of 
a discrete parity symmetry. Thus, the neutral components of the extra doublet are stable and 
can be considered as probable candidate of Cold Dark Matter. We have made calculations for 
three different types of Dark Matter experiments, namely, 76 Ge (like GENIUS), DAMA (Nal) 
and XENON ( 131 Xe). Also demonstrated the annual variation of Dark Matter detection in 
case of all three detectors considered. 
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A 'Heavy Higgs Model' has been proposed [1] within the framework of SU{2) L x U(1)y 
symmetry through the inclusion of another doublet and a discrete parity symmetry. This 
model is motivated to give an alternate way to control the quadratic divergence in the Higgs 
sector and to solve the naturalness problem. Again, results of electroweak precision tests 
are in favour of light Standard Model (SM) Higgs boson , m h < 186 GeV at 95% c.l [2] 
with a central value too low from the lower bound (114 GeV) obtained from direct search 
experiments whereas, the present model is based on adopting the few per cent chance that 
there is no Higgs boson mass below 200 GeV. The basic feature of the model is as follows. 

The additional doublet neither develops VEV nor it has any coupling with matter, how- 
ever, the doublet reprsentation necessarily admits weak interaction and also scalar self in- 
teraction. Due to absence of any coupling with matter of the extra doublet its neutral 
components are stable, and, therefore, those Lightest Stable Higgs (LSH) may be considered 
as probable candidate of cold dark matter. The zero VEV of the extra doublet is assured 
through the choice of coefficients of the scalar potential of the model and the discrete parity 
symmetry proposed in the model prohibits the extra doublet to couple with fermions. 

In the present work, we investigate Direct detection rates of those Dark Matter particles 
in view of ongoing and future Dark Matter experiments. 

Let us first discuss relevant parts of the model. The model contains two doublet scalars Hi 
and H 2 . The H\ scalar is the usual standard model doublet and the additional H 2 scalar 
is fermiophobic. The discrete parity symmetry imposed to achieve decoupling of H 2 field is 
given by 

H 2 - -H 2 (1) 
while keeping all other fields invariant. The VEV of the Higgs fields are 

H 1 = (0,v), # 2 = (0,0) (2) 

and expanding around their minima we get 




The scalar potential of the model is given by 



V = ^|// 1 | 2 + ^|// 2 | 2 + A 1 | J H" 1 | 4 + A 1 |// 2 | 4 + A 3 |//i| 2 |// 2 | 2 

+\,\H\H 2 \ 2 + ^ [{H\H 2 f + h.c] . (4) 
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The mass of the usual Higgs boson h in the present model is set around m h ~ 400 — 600 
GeV and the three additional scalars as 

m] = n\ + A/w 2 , 1 = H,S,A 

= A3 

A5 = A3 + A4 + A5 

Aa = A 3 + A 4 - A 5 (5) 

and from the minimization of the scalar potential one obtains conditions for which the 
potential is stable, as 

Ai, 2 > and A 3 , X L > -2(AiA 2 ) 1/2 (6) 

where 

A L = A3 + A4- |A 5 | (7) 

The Higgs fields S, A are stable as long as the parity symmetry (1) is unbroken and those 
stable neutral scalars appear to be a candidate for Cold Dark Matter in the universe. In 
order to see the prospect of direct detection of such a Dark Matter candidate, we consider 
two processes of scalar-nucleon interaction i) via Z boson exchange and ii) via Higgs scalar 
h exchange. 

The amplitude comes out due to Z boson exchange is too large (about 8-9 orders) com- 
pared to the existing CDMS collaboration experimental limits [3] on spin-independent wimp- 
nucleon interactions from the two-tower run of the cryogenic dark matter search. Hence, to 
neglect this contribution, a mass splitting between the two scalars S and A are considered 
which is greater than the kinetic energy of the dark matter in our galactic halo, so that the 
process due to Z exchange is kinematically forbidden. We consider the other process due to 
h exchange at the tree level for our analysis in the present work. 

In Ref. [1], the possible range of mass, m L for the LSH, to account for the Dark Matter 
content of the universe, is discussed. For m L <; m w (mass of W boson), the Dark Matter 
density Odm falls much below the acceptable value due to dark matter annihillation to W 
bosons [LL — ► WW, where L is generic representation of S and A fields). Therefore, 
massive LSHs can contribute to only a fraction of Dark Matter content. While LSH mass 
m,L & rn w , the annihillation to W bosons are kinematically disfavoured and in particular 
m L ~ 60 — 70 GeV can account for the Dark Matter content of the universe with Am, 
the mass difference between the LSH and NLSH (next lightest stable Higgs), ~ 8 GeV. 
Therefore, we consider three LSH masses, in the range 60 - 70 GeV in the present work. We 
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also consider, as an example, a value of in the ttil » my/ regime, which as mentioned 
before, can only account for just a fraction of the total Dark Matter. 



The scalar-nucleon cross-section due to the later case comes out as 




m 2 . ( \l 



) 



2 



a(LN -> LN) 



4tt \mLm\ 




(8) 



where m r is the reduced nucleon mass, Xl is the combination of scalar couplings defined in 
eqn.(4), m h is the mass of the Standard Model like Higgs scalar and f N is the usual nuclear 
matrix element given by 



and for the present analysis we set fx ~ 0.3. 

In the present work, we calculate the possible direct detection rates for such Lightest Stable 
Particles (or LSH) Dark Matter candidates discussed above, in the experiments like GENIUS 
(target material 76 Ge) [4, 5], DAMA (target material Nal) [6, 7, 8] and XENON (target 
material 131 Xe) [9, 10]. 

The direct detection of Dark Matter with a terrestrial detector uses the elastic scattering 
of Dark Matter candidate off the detector nuclei. As this cross-section is very small, the 
energy deposited by a Dark Matter candidate of mass in the range 1 GeV to 1 TeV on a 
detector nucleus is not generally more than 100 keV. Hence to perform this difficult task of 
Dark Matter detection a very low threshold detector condition is required. 

Differential detection rate of Dark Matter per unit detector mass can be written as 



where Nt denotes the number of target nuclei per unit mass of the detector, $ - the Dark 
Matter flux, v - the Dark Matter velocity in the reference frame of earth with f(v) - its 
distribution. The integration is over all possible kinematic configurations in the scattering 
process. In the above, |q| is the momentum transferred to the nucleus in Dark Matter- nucleus 
scattering. Nuclear recoil energy E R is 



< N\Em q qq\N >= f N m N < N\N > 



(9) 




(10) 



E R = \q\ 2 /2m- 



nuc 



m% ed v 2 (l — COS 9) j TTlj 

m L m nuc 



nuc 



(11) 

(12) 



m red = 



m L + m DUC 
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where 9 is the scattering angle in Dark Matter-nucleus centre of momentum frame, m nuc is 
the nuclear mass and m L is the mass of the Dark Matter. 

Now expressing $ in terms of local Dark Matter density p x , velocity v and mass mi and 
writing |q| 2 in terms of nuclear recoil energy En with noting that Nt = l/m nuc , Eq. (10) 
takes the form 



dR 



p x da 



dE 



R 



111 1, d (J ./ r„ .. 



2 

,1/2 



vf(v)dv, 



2m 2 ed 



(13) 



Following Ref. [11] the Dark Matter- nucleus differential cross-section for the scalar in- 
teraction can be written as 

If = i& 2 ™ • (14) 

In the above cr sca i ar is Dark Matter-nucleus scalar cross-section and F{Er) is nuclear form 
factor given by [12, 13] 



F(E R ) 



cxp 



(15) 



qRi 

Rl = (r 2 -5s 2 ) 1/2 
r = 1.2A 1/3 

where thickness parameter of the nuclear surface is given by s ~ 1 fm, A is the mass number 
of the nucleus and j\{qR\) is the spherical Bessel function of index 1. 

The distribution /(%a) of Dark Matter velocity t> ga i with respect to Galactic rest frame, 
is considered to be of Maxwellian form. The velocity v (and f(v)) with respect to earth rest 
frame can then be obtained by making the transformation 

V = V ga i - V e (16) 

where ■% is the velocity of earth with respect to Galactic rest frame and is given by 

/ 2vr(t - t ) 



Vq + IVb cos 7 COS 



17) 



In Eq. (17), T = 1 year, the time period of earth motion around the sun, to — 2 nd June, t> or b 
is earth orbital speed and 7 ~ 60° is the angle subtended by earth orbital plane at Galactic 
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plane. The speed of solar system v Q in the Galactic rest frame is given by, 



Vq = V + V 



pec 



(18) 



where v is the circular velocity of the Local System at the position of Solar System and v pec 
is speed of Solar System with respect to the Local System. The latter is also called peculiar 
velocity and its value is 12 km/sec. Although the physical range of vq is given by [14, 15] 
170 km/sec < v < 270 km/sec (90 % C.L.), in the present work we consider the central 
value of Vq = 220 km/sec. Eq. (17) gives rise to annual modulation of Dark Matter signal 
reported by DAMA/Nal experiment [6, 7, 8]. 

Defining a dimensionless quantity T(E R ) as, 



nE R ) = r 

and noting that T(E R ) can be expressed as [11] 



/(«) 



dv 



(19) 
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R 



-Vq 



erf 



+ 



we obtain from Eqs. (13) and (14) 

dR ^scalarPx 



dE R 



F 2 {E 



— erf 



(20) 



4w e m L m r 2 cd 



R 



erf 



+ 



(21) 



The total local Dark Matter density p x is taken to be 0.3 GeV/cm 3 . The above expression 
for differential rate is for a monoatomic detector like Ge but it can be easily extended for a 
diatomic detector like Nal as well. 

The measured response of the detector by the scattering of Dark Matter off detector 
nucleus is in fact a fraction of the actual recoil energy. Thus, the actual recoil energy E R is 
quenched by a factor qnx (different for different nucleus X) and we should express differential 
rate in Eq. (21) in terms of E = qn x E R . For 76 Ge, qn Ge = 0.25 [16], for f a , qn Na = 0.3 [17], 
for 127 I, q ni = 0.09 [17] and for 131 Xe, qn Xe = 0.8 [16]. 

Thus the differential rate in terms of the observed recoil energy E for a monoatomic 
detector like Ge detector can be expressed as 



AR 

~Ke 



(E) 



(E+AE)/ q n Ge dR G dE R 

E/ g n Gc dE R AE 



(22) 
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and for a diatomic detector like Nal, the above expression takes the form 

AR , r{E+AE)/ q n^ dR ^ dE R 

AE {E) ~ aNa A MNa 1e^ {Er) AE 

r(E+AE)/ qni d Rl d ER 



+ ai 



JE /qm dE- R {ER) AE (23) 



where ciN a and a\ are the mass fractions of Na and I respectively in a Nal detector and are 
given by (see Table 2) 

aNa = = 0.153 ai = = 0.847 

™ Na + m m Na + mi 

The differential detection rates AR/ AE (/kg/day/keV) in case of LSH Dark Matter for 
different values of observed recoil energies are calculated using Eqs (10 - 22) for monoatomic 
detectors like Ge and Xe and using Eqs (10 - 23) for diatomic detector like Nal with AE = 1 
keV. For calculation of the differential rates, we put t — t in Eq. (17) for all three types of 
detectors considered here. 

In Fig. 1 we plot LSH Dark Matter detection rate for 76 Ge detector such as in GENIUS 
experiment at Gran Sasso. For demonstrative purpose we calculate the rates for three 
different LSH masses in m L < raw regime, namely m L = 60 GeV, 65 GeV and 70 GeV. 
For all the calculations the value of the coupling |A| is fixed at 0.5. It appears from Fig. 1 
that the differential rate decreases with the increase of LSH mass. This can be understood 
from Eq. (21) where the LSH mass mi appears at the denominator. We wish to make a 
comment here that while trying to plot this variation in logscale, it reveals that although 
the rates are more for lower mass of LSH Dark Matter and decreases with the increase of 
mass (as is evident form Fig. 1), the situation becomes reversed for high recoil energies as 
the rates show an oscillatory behaviour which becomes prominent at higher recoil energies. 
This phenomenon is due to the oscillatory nature of the Bessel function used in Eqs. (14) 
and (15). But at those high recoil energies where the oscillatory nature becomes prominent, 
the yields are virtually nil. Hence we do not plot them in linear scale instead. 

The variation of rates with detector recoil energies for diatomic Nal detector (used in 
DAMA experiment), for three different LSH masses mentioned above, is shown in Fig. 2. 
The coupling |A| is fixed at the same value of 0.5. 

Similar results for Xenon ( 131 Xe) detector are plotted in Fig. 3. 

One very positive signature of Dark Matter detection by direct detection method is 
the periodic annual variation of the detected Dark Matter. This periodicity arises due 
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to the periodic motion of the earth around the sun by which the directionality of earth's 
motion changes continually throughout the year. Due to this, the amount of Dark Matter 
encountered by the earth varies annually which can be understood from Eq. (17). This 
is the annual modulation of the detected Dark Matter and detecting this annual variation 
serves as a confirmatory test for Dark Matter detection. In order to see the possible annual 
variation of the detected LSH Dark Matter in three types of detectors discussed here, we 
have calculated the total events per day for a whole year for each of these three types of 
detectors. Thus the value of t in Eq. (17) varies from 1 to 365 while to — 153 (2 nd June). 
For this purpose we have chosen a LSH mass of 65 GeV and the value of coupling to be 
|A| = 0.5. The results are plotted in Figs. 4, 5 and 6 for 76 Ge, Nal and 131 Xe detectors 
respectively. The plots clearly show the sinusoidal behaviour of daily yield over a year. It 
peaks in June and is minimum in December as expected. 

For comparison, we calculate the rate and annual variation for one sample case for » 
raw regime and as discussed earlier this LSH cannot account for the total Dark Matter. Hence 
the local Dark Matter density p x for this case is taken to be 0.003 GeV/cm 3 rather than the 
usual total local Dark Matter density of 0.3 GeV/cm 3 . The coupling |A| for this case is also 
kept at 0.5. Both the rate (Fig. 7) and the annual variation (Fig. 8) for this LSH clearly 
orders of magnitude lower than the cases considered in mj,^ mw regime (Figs 1-6). 

In summary, we investigate direct detection rates for a possible Dark Matter candidate 
in a Heavy Higgs model. The model contains an additional Higgs which neither couples 
with matter nor it has developed any VEV. This can be achieved through the choice of 
model parameters and discrete parity symmetry. Thus, the neutral part of the extra doublet 
becomes stable and can be a possible candidate for Cold Dark Matter. The LSH with mass 
around 60 - 70 GeV (in the regime mi ^ mw) can explain the Dark Matter content of the 
universe. On the other hand if mi > mw, the LSH candidate may explain only a fraction 
of Dark Matter. We calculated direct detection rates of those Dark Matter candidates in 
the context of three experiments namely 76 Ge (as in GENIUS), DAMA and XENON. We 
have also shown the annual variation of the Dark Matter detection due to periodic motion 
of earth, for all the three experiments considered here. For comparison, the results for one 
case with m L = 300 GeV (>> m w ) is also shown. 
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Figure Captions 



Fig. 1 Differential detection rates AR/ AE vs recoil energy E for three different values of 
LSH mass (m L ) namely 60 GeV, 65 GeV and 70 GeV. The topmost plot corresponds to 
mi = 60 GeV and the lowermost plot corresponds to m L = 70 GeV and the plot in between 
is for rriL = 65 GeV. The coupling constant is kept fixed at 0.5 (see text). 

Fig. 2 Same as Fig. 1 but for Nal. 

Fig. 3 Same as Fig. 1 but for 131 Xe. 

Fig. 4 Annual variation of LSH Dark Matter direct detection signal (events/kg/day vs each 
day in a year) for 76 Ge detector with tjil = 65 and |A| = 0.5. 

Fig. 5 Same as Fig. 4 but for Nal. 

Fig. 6 Same as Fig. 4 but for 131 Xe. 

Fig. 7 Differential detection rates AR/ AE vs recoil energy E for LSH mass m L = 300 GeV 
for 76 Ge. 

Fig. 8 Annual variation for LSH mass m L = 300 at 76 Ge detector. 
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